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TIME VARIANT FILTER IMPLEMENTATION 
Field of the Invention 

[0001] The present invention relates to filtering, and 

in particular to implementing a time variant filter 
cascaded with a time invariant filter. 

Backgroxxnd of the Invention 

[0002] In many communication applications, one is 

required to realize a time variant mult i -rate 
interpolation filter followed by a time invariant filter, 
or a time invariant filter followed by a time variant 
multi-rate decimation filter. For these two filter 
configurations, the direct realization of the time 
invariant filter is computationally intensive because it 
operates at a higher sampling rate as compared to the 
sample rate of the input of the interpolation filter or 
the output of the decimation filter. In order to reduce 
the computational cost of the fixed filter, the time 
invariant filter can be combined with the time variant 
multi-rate interpolation/decimation filter in such a 
manner as to have the combined filter effectively operate 
at the lower sampling rate. 

[0003] Such filtering techniques find practical 
application in the area of multiple access interference 

(MAI) mitigation. In direct sequence code-division 
multiple access (DS-CDMA) systems, MAI is a factor that 
contributes to a limitation in system capacity and 
performance. In an attempt to reduce the effect of this 
factor, one can employ some form of multi-user detection 

(MUD) algorithm. The basis of a typical MUD detector is 
the application of information known about other users to 
improve detection of each individual user. 

[0004] Of particular interest is the class of MUD 

detectors known as subtractive interference cancellation 
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detectors. The fundamental principle behind these 
detectors is that an estimate is made of each individual 
user's contribution to the total MAI and then subtracted 
out from the received signal such that the MAI affecting 
each individual user is reduced. Typically, the detection 
process is carried out in an iterative manner, where the 
data decisions of the previous iteration are used as the 
basis for the next iteration's MAI estimates. 
[0005] For a chip-rate interference cancellation 
detector, the process of calculating each individual's 
contribution to the total MAI involves "regenerating" 
chip-rate samples from symbol rate data estimates. As an 
example, for the time division duplex (TDD) mode of a 
Universal Mobile Telecommunications System terrestrial 
radio access (UTRA) communication system, regeneration 
for a given user involves: spreading by the user's 
corresponding real channelization code, scrambling with a 
cell-specific complex code, and filtering with the user's 
channel estimate, which typically requires 2 8 to 64 
complex taps. Note that the spreading codes generally 
vary in length from 1 to 16. However, the cell-specific 
scrambling code usually has a fixed length of 16. 
Therefore, if spreading and scrambling are combined into 
a single function, the result is effectively a time 
variant multi-rate interpolation filter. In this way, the 
regeneration process can be considered as a time variant 
multi-rate interpolation filter cascaded with a time 
invariant filter. 

[0006] Given that regeneration is performed on a per 

cancellation iteration basis for each user, the 
computational cost of regeneration is substantial. As 
such, there is a need to reduce the computational cost of 
regeneration without degrading the accuracy of the 
regenerated results. Conventionally, regeneration is 
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accomplished by first performing channelization in the 
traditional manner, such as by effecting a single, time 
variant complex multiplication at the chip rate. Finite 
impulse response (FIR) filters, whose complex 
coefficients are given by each user's corresponding 
channel estimate, follow. Since these FIR filters are 
operating at the chip-rate, the computational 
requirements are very high. 

[0007] As an example, consider a TDD mode UTRA system 
with eight active users and channel estimates with 57 
complex taps. The total computational cost of 
regeneration would be 4 64 complex calculations per chip- 
rate period, per cancellation iteration. Assuming that 
interference cancellation uses four cancellation 
iterations per TDD slot period (equivalent to 2560 chip- 
rate periods or 667 microseconds) , this corresponds to a 
computational requirement of 7.2 billion complex 
operations per second. An alternative technique that has 
been proposed involves the cancellation of interference 
at the symbol -rate as opposed to the chip-rate. This 
approach eliminates the need for chip-rate regeneration, 
but also may result in a reduction in the accuracy of the 
MAI estimates. 

[0008] Accordingly, there is a need for efficient 
realization of a time variant multi-rate interpolation, 
or decimation FIR filter, cascaded with a time invariant 
FIR filter. 

Summary of the Invention 

[0009] The present invention provides efficient 

realization of a time variant multi-rate FIR filter 
cascaded with a time invariant FIR filter. Preferably, 
the time invariant FIR filter response provides filtering 
based on channel estimates and the time variant FIR 
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filter provides filtering to effectively implement 
spreading and scrambling with appropriate coding. The 
realization of the combined filters is accomplished by 
transforming the pair of filters into a single "combined" 
polyphase filter with coefficients that are time variant 
in a specific manner. The invention provides for pre- 
calculating the coefficients of the combined filter; a 
memory storage and addressing scheme for the resulting 
coefficients that ensures that the combined filter 
behaves like the original pair of filters; and overall 
polyphase decomposition of the resulting filter to 
significantly reduce the computational complexity as 
compared to the requirements of the original pair of 
filters . 

[0010] Those skilled in the art will appreciate the 

scope of the present invention and realize additional 
aspects thereof after reading the following detailed 
description of the preferred embodiments in association 
with the accompanying drawing figures . 

Brief Description of the Drawing Figures 

[0011] The accompanying drawing figures incorporated 

in and forming a part of this specification illustrate 
several aspects of the invention, and together with the 
description serve to explain the principles of the 
invention. 

[0012] Figure 1 is a schematic representation of a 

base station transceiver configured according to one 
embodiment of the present invention. 

[0013] Figure 2 is a schematic representation of 

multiple access interference cancellation circuitry 
according to one embodiment of the present invention. 

[0014] Figure 3 is a schematic representation of a 

time variant filter preceding a time invariant filter. 
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[0015] Figure 4 is an illustrative memory architecture 
of combined time variant and time invariant filters for 
the architecture illustrated in Figure 3 . 

[0016] Figure 5 is a schematic representation of a 
time variant filter following a time invariant filter. 

[0017] Figure 6 is an illustrative memory architecture 
of combined time variant and time invariant filters for 
the architecture illustrated in Figure 5 . 

[0018] Figure 7 is a representation of a time variant 

multi-rate interpolation filter. 

[0019] Figure 8 is a representation of a time variant 
multi-rate decimation filter. 

[0020] Figure 9 is a representation of a fourth-order 
combined filter. 

[0021] Figure 10 is an illustrative memory array for a 

fourth-order filter. 

[0022] Figure 11 is a table of example filter outputs 
for a fourth-order combined filter. 

[0023] Figure 12 is a representation of a fourth-order 

interpolation by three filter. 

[0024] Figure 13 is a table of example filter outputs 

for a fourth-order interpolation by three filter. 
[0025] Figure 14 is a representation of a fourth-order 
polyphase interpolation by three filter. 

[0026] Figure 15 is a table of example filter outputs 

for a fourth-order polyphase interpolation by three 
filter. 

[0027] Figure 16 is a representation of a regeneration 

process including spreading and scrambling functions 
implemented in hardware. 

[0028] Figure 17 is a representation of a regeneration 
process including spreading and scrambling functions 
implemented using a time variant filter. 
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Detailed Description of the Preferred Embodiments 
[0029] The embodiments set forth below represent the 

necessary information to enable those skilled in the art 
to practice the invention and illustrate the best mode of 
practicing the invention. Upon reading the following 
description in light of the accompanying drawing figures, 
those skilled in the art will understand the concepts of 
the invention and will recognize applications of these 
concepts not particularly addressed herein. It should be 
understood that these concepts and applications fall 
within the scope of the disclosure and the accompanying 
claims . 

[0030] The present invention relates to the efficient 
realization of a time variant multi-rate finite impulse 
response (FIR) filter cascaded with a time invariant 
multi-rate FIR filter, the overall result of which could 
implement interpolation, decimation, or a combination 
thereof. This is accomplished by transforming the pair of 
filters into a single "combined" polyphase filter with 
coefficients that are time variant in a specific manner. 
The invention provides for pre-calculating the 
coefficients of the combined filter; a memory storage and 
addressing scheme for the resulting coefficients that 
ensures that the combined filter behaves like the 
original pair of filters; and overall polyphase 
decomposition of the resulting filter to significantly 
reduce the computational complexity as compared to the 
requirements of the original pair of filters. In the 
preferred embodiment, the invention finds practical 
application in the efficient calculation and estimation 
of multiple access interference (MA.I) in an interference 
cancellation detector for direct -sequence code-division 
multiple access (DS-CDMA) communication systems. Those 
skilled in the art will recognize other applications of 
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the invention. The following description initially 
defines the operating environment of the present 
invention . 

[0031] The present invention may be incorporated in a 
base station transceiver 20. The basic architecture of a 
base station transceiver 20 is represented in Figure 1 
and may include a receiver front end 22, a radio 
frequency transmitter section 24, an antenna 26, a 
duplexer or switch 28, a baseband processor 30, a control 
system 32, and a frequency synthesizer 34. The receiver 
front end 22 receives information bearing radio frequency 
signals from one or more remote transmitters provided by 
a terminal, such as a mobile telephone, wireless personal 
digital assistant, or like communication device. A low 
noise amplifier 38 amplifies the signal. A filter 
circuit 40 minimizes broadband interference in the 
received signal, while downconversion and digitization 
circuitry 42 downconverts the filtered, received signal 
to an intermediate or baseband frequency signal, which is 
then digitized into one or more digital streams. The 
receiver front end 22 typically uses one or more mixing 
frequencies generated by the frequency synthesizer 34. 
[0032] The baseband processor 30 processes the 

digitized received signal to extract the information or 
data bits conveyed in the received signal . This 
processing typically comprises demodulation, decoding, 
error correction, and inference cancellation operations. 
As such, the baseband processor 3 0 is generally 
implemented in one or more digital signal processors 
(DSPs) , application specific intergrated circuits 
(ASICs) , or field programmable gate arrays (FPGAs) . 
Further detail regarding the operation of the baseband 
processor 3 0 is provided below. On the transmit side, 
the baseband processor 30 receives digitized data, which 
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may represent voice, data, or control information, from 
the control system 32, which it encodes for transmission. 
The encoded data is output to the transmitter 24, where 
it is used by a modulator 44 to modulate a carrier signal 
that is at a desired transmit frequency. Power amplifier 
circuitry 46 amplifies the modulated carrier signal to a 
level appropriate for transmission, and delivers the 
modulated carrier signal to antenna 2 6 through a matching 
network 4 8 . 

[0033] In a typical parallel interference cancellation 

(PIC) system addressing MAI, the input signal is composed 
of multiple (M) user signals separated by orthogonal 
codes. For a specific time period, or slot, the input 
signal is separated into M individual signals 
corresponding to each of the multiple users and 
demodulated to recover K symbols for each of the M 
individual signals. In a "regeneration" process, the 
respective K symbols for the M users are remodulated 
using the appropriate spreading and scrambling coding and 
filtered using channel estimates to create M regenerated 
signals. For each user, the regenerated signals for the 
other M-1 users are subtracted from the original input 
signal to form M new individual signals. This last step 
is repeated for each user to create M individual signals 
for each of the M users. 

[0034] These new individual signals are individually 
demodulated and reprocessed, ideally with higher 
reliability than in the previous step, for as many 
iterations as desired. The output of the last iteration 
is used to form the final symbol decisions from the each 
of the M user signals. A survey of interference 
cancellation techniques for CDMA is given in S. Moshavi, 
"Multi-User Detection for DS-CDMA Communications," IEEE 
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Communications Magazine, October 1996, pp. 124-136, which 
is incorporated herein by reference. 

[0035] The architecture representative of PIC detector 
circuitry for one embodiment of the present invention is 
illustrated in Figure 2. This particular variant of PIC 
is referred to as "chip-level" PIC due to the fact that 
interference is cancelled at the chip-rate. While chip- 
rate cancellation is computationally intensive, it leads 
to a simple data path-oriented structure that remains 
relatively constant regardless of the given number of 
users and spreading factor combinations. The baseband 
processor 3 0 will preferably include the circuitry and 
software to implement the following structure or 
functionality. Those skilled in the art will recognize 
the numerous alternatives to the specifically disclosed 
structure and function. In general, the following 
functions are provided in hardware, software, or a 
combination thereof: midamble extraction 62, channel 
estimation 64, parallel interference cancellation 66, 
rake demodulation 68, symbol decision 70, symbol 
processing 72, and signal regeneration 74. Signal 
regeneration 74 is significantly enhanced by the present 
invention . 

[0036] For the disclosed embodiment, it is assumed 

that a single chip-rate PIC application specific 
integrated circuit (ASIC) is used to process all users 
for a single antenna. Moreover, it is assumed that the 
over-sampled, received signal is decimated to the chip- 
rate so that the input signal 60 is a composite signal 
represented by chip-rate complex samples. It is also 
assumed that channel estimation and other required pre- 
calculations occur during the reception of the second 
data field for a given time slot, and the processing of 
all data symbols for a given time slot occurs during the 
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reception of data for the next time slot. Those skilled 
in the art will recognize equivalent implementation as 

those specifically described herein. 

[0037] Initially, the midamble extraction function 62 

extracts the midamble from the received input signal 60 
for a given slot, and a channel estimate is obtained for 
each user from the channel estimation function 64 in 
traditional fashion. The channel estimates may be used 
for chip-rate finger allocation in the constituent rake 
demodulator 68 and used to directly or indirectly define 
the filter coefficients for the regeneration function 74. 
Channel estimation may be accomplished using a "fast" 
cyclic correlator and additional processing to enhance 
the accuracy of the result . 

[0038] Upon reception, the PIC function 66 is bypassed 
and rake demodulation 68 separates the composite input 
signal 6 6 into M individual signals corresponding to the ■ 
respective signals for each user. To determine the 
individual signals for each user during subsequent 
iterations, the PIC function 66 subtracts the regenerated 
signals corresponding to all other individual signals 
from the composite input signal 66 for the given time 
slot. The regenerated signals for each user are fed back 
via the regeneration and symbol processing functions 74, 
72 . 

[0039] As noted, rake demodulation 68 is used to 
demodulate each individual user' s signal from the 
composite input signal 60. Demodulation may include 
despreading each individual signal with a unique 
channelization code and descrambling all individual 
signals using the base station's descrambling code. This 
operation is preferably performed at the chip-rate, 
producing a complex, symbol-rate output. Symbol estimates 
for each user are determined by applying the symbol 
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decision function 70 to the outputs of the rake 
demodulator 68 in traditional fashion. The symbol 
decision function 70 may be implemented in hardware using 
symbol detection circuitry. 

[0040] Symbol processing 72 may be performed on the 
symbol estimates in order to enhance cancellation 
performance. A preferred method of processing symbols is 

provided in U.S. patent application number / , 

filed , entitled SYMBOL-DIRECTED WEIGHTING IN 

PARALLEL INTERFERENCE CANCELLATION. The symbol 
processing 72 in the preferred method effectively 
provides a weighting for the latest symbol decision based 
on the consistency of previous estimates for the given 
symbol in earlier iterations. The symbol processing 72 
provides individually weighted symbol estimates for each 
symbol based on the weighting and symbol estimate. 
[0041] The regeneration function 74 essentially 

recreates estimates of the individual signals from the 
detected symbols, or weighted symbol estimates if such 
processing is used, for each of the estimated users. 
Accordingly, regeneration will typically involve 
processing to effectively: 

-spread the symbols weighted with the user's 
specific channelization or spreading code, 

-scramble the spread signal with a base station's 
complex code, and 

-filter the spread and scrambled signal with the 
user's channel estimates. 

[0042] The output of the regeneration function 74 

provides estimates for each of the individual 
contributions to the composite received signal received 
for each user. As described in detail below, the present 
invention reduces the computational complexity of the 
regeneration operation without introducing any error. 
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Each of the regenerated signals is individually sent to 
the PIC function 66 to effect the iterative MAI 

cancellation . 

[0043] For the regeneration function 74, the present 
invention efficiently cascades a time variant FIR filter 
and a time invariant FIR filter and effects a polyphase 
decomposition of the result. The polyphase decomposition 
leads to a reduction in complexity because the resulting 
polyphase filter can effectively operate at a lower 
sampling rate. 

[0044] In essence, the invention combines a time 
variant filter and time invariant filter into one time 
variant polyphase filter that provides the same transfer 
function as separately implemented filters. The present 
invention is not limited to regeneration in a PIC 
environment and can be used for filter applications that 
match the general configurations of Figures 3 and 5. 
Further, the ordering of the filters being combined may 
be changed and the filters may be implemented in 
interpolation and decimation embodiments. In general, 
the invention includes: 

- a pre-calculation technique for generating all 
possible coefficients of the combined time variant 
polyphase filter; 

- a memory storage and addressing scheme for the 
resulting filter coefficients such that the appropriate 
transfer function is maintained; and 

- overall polyphase decomposition of the resulting 
filter to reduce the computational complexity. 

[0045] The particular technique for combing the two 

FIR filters and the associated pre-calculation and 
storage of the resulting coefficient differs depending on 
whether the time variant filter is effectively 
implemented before or after the time invariant filter. 
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Figure 3 illustrates a time variant filter hv(n) preceding 
a time invariant filter hf (n) . Let the input-output 
relationship of the time variant filter hv(n) be 

Eq. 1 

4=0 

and the input -output relationship of the time invariant 
filter hf(n) be 

Eq. 2 

1=0 

where and Nf are the orders of the time variant and 
time invariant FIR filters, respectively. From Eq. 1 and 
Eq. 2, the input -output relationship of the combined 
filter in Figure 3 can be reduced to 



Eq. 3 



yin)= 2] 



x(n - m) _ 



Then, one can simplify Eq. 3 in accordance with 
Eq. 4 



which has time variant coefficients c^{n) given by 
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Eq. 5 

mm(m,iV„) 
*=max(0,m-;v^) 

[0046] In this way, one obtains a single combined time 

variant filter. Based on Eq. 5, the combined filter 
coefficients c„(n)(m e (0,//^ +iV^)) are not only a function of 
the different coefficients of the time variant filter, 
but also a function of these coefficients at different 
points in time. Therefore, even if all coefficients of 
the time variant filter change at the same time, future 
values of Cn,{n) will be affected by this change. As such, 
it takes time for the coefficient changes to "ripple" 
through and affect all of the combined filter 
coefficients. If certain assumptions are made about the 
characteristics of the time variant filter coefficients 
akin) , then a simplified architecture for the resulting 
combined filter may be created. These assumptions are as 
follows: a) there are R sets of unique values for the 
variable coefficients akin), and b) the values of akin) 
are changed at the same time every P sample periods. If 
these assumptions hold, the resulting unique values of 
Cra(n) will be finite and will change in a periodic 
fashion. The period of this variation is given by 

Eq. 6 

'^VARIATION = P • 

Note that for a given m, the number of unique values for 
Cmin) will be less than or equal to Tvariation, depending on 
the values of m, R, P, N^, and Nf. 
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[0047] Using these results, one can obtain a 
simplified memory architecture for the hardware 

realization of the combined filter. Coefficient values 
are stored in a memory array with Tvariation rows and + Nf 
+ 1 columns. Each combined filter coefficient 
c^(n)(m e(0,N^ + Nf)) is allocated a corresponding column in 
this memory array. The values of c^{n) are pre-calculated 
in accordance with Eq . 5, taking into account the 
possible values for a^di) , and stored in the corresponding 
column. Note that in order for the set of coefficient 
values for a given sample period to be contained in the 
same row, the values of Cmin) must be calculated and 
stored in a specific manner. An example of an algorithm 
for combined coefficient calculation follows. 



Algorithm 1 

//Initialize memory array 
for m = 0 to Nv + Nf 

for J = 0 to Tvariation - 1 
c„,(j) = 0 

end 

end 

//Perform combined coefficient calculation 
for m = 0 to + Nf 

for k = max (0, in - Nf) to min {m, N^) 

index = mod {m-k, Tvariation) 

for i = 0 to i? - 1 

for J = 0 to P - 1 

Cmi index) = Cmi index) + ak{i*P) *bin-k 
index = mod ( index + 1 , Tvariation) 

end 

end 

end 

end 

[0048] Using such an algorithm allows one to use a 
common memory index in order to select the coefficient 
values of the combined filter for a given sample period. 
During normal operation of the combined filter, the 
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memory index is initialized to row 0 and incremented by- 
one row (modulo Tvariation) every sample period. A graphical 
depiction of the memory architecture of the resulting 
filter is given in Figure 4 . Once the combined filter 
has been obtained, one can then perform a conventional 
polyphase decomposition of the result in order to realize 
the computational savings. 

[0049] Figure 5 illustrates a time variant filter hv(n) 
following a time invariant filter hf (n) . As described 
above, let the transfer function of the time variant 
filter hv(n) be represented by Eq. 1 and the transfer 
function of the time invariant filter hf(n) be represented 
by Eq. 2. From Eq. 1 and Eq, 2, the transfer function of 
the combined filter in Figure 5 can be represented as 



Eq. 7 




k=0 1=0 



J 



Then, one can simplify Eq, 7 in accordance with 



Eq. 8 



m=0 



which has time variant coefficients Cmiri) given by 



Eq. 9 



k=max.(fi,m-Nf) 
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[0050] In this way, one obtains a single combined time 

variant filter. From Eq. 9, it can be seen that the 
combined filter coefficients c^{n){m^{0,N^+Nj-y) are only a 
function of the different coefficients of the time 
variant filter at the same time instant. Therefore, if 
all coefficients of the time variant filter change at the 
same time, all values of Cmin) will be immediately 
affected by this change. Accordingly, changes in 
coefficients occur immediately and do not take time to 
"ripple" through the memory to affect all of the combined 
filter coefficients. 

[0051] As before, if certain assumptions are made 
about the characteristics of the time variant filter 
coefficients 3^(21), then a simplified architecture for the 
resulting combined filter may be created. These 
assumptions are as follows: a) there are R sets of unique 
values for the variable coefficients akin), and b) the 
values of akin) are changed at the same time every P 
sample periods. If these assumptions hold, the resulting 
unique values of c^in) will be finite and will change in a 
periodic fashion with a period given by Eq. 6. 

[0052] Given that changes to a^in) immediately affect 
all coefficient values c^in) , the number of unique values 
for each combined coefficient is exactly R. Using these 
results, one can obtain a simplified memory architecture 
for the hardware realization of the combined filter. 
Coefficient values are stored in a memory array with J? 
rows and + Nf + 1 columns. Each combined filter 
coefficient c„(n){m g (0,N^ + Ny)) is allocated a 
corresponding column in this memory array (i? memory 
locations) . The values of coefficients c^(n) are pre- 
calculated in accordance with Eq. 9, taking into account 
the possible values for a^Cn) , and stored in the 
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corresponding column. An example of an algorithm for 
combined coefficient calculation for this case follows. 

Algorithm 2 

//Initialize memory array 
for m = 0 to + Nf 

for j = 0 to i? - 1 

C:r,ij) = 0 

end 

end 

//Perform combined coefficient calculation 
for m = 0 to + Nf 

for k = maKiO,m - Nf) to minim, N^) 
for i 0 to i? - 1 

Cmii) = Cm(i) + ak{i*P) * b^-k 

end 

end 

end 

[0053] In this way, a common memory index can be used 
to select the coefficient values of the combined filter 
for a given sample period. During normal operation of 
the combined filter, this memory index is initialized to 
row 0 and incremented by one row (modulo R] every P 
sample periods. A graphical depiction of the memory 
architecture for the resulting filter is given in Figure 
6. Once the combined filter has been obtained, one can 
then perform a conventional polyphase decomposition of 
the result in order to realize the computational savings. 
[0054] Polyphase decomposition corresponds to a 
partitioning of the overall combined time variant filter 
into a set of "sub-filters" arranged in a specific 
manner. The coefficient values for these sub- filters are 
directly obtained from the current row, as indicated by 
the memory index, in the pre-calculated coefficient 
memory array. The resulting filter structures are given 
in Figure 7 and Figure 8, for the case of a time variant 
interpolation filter and a time variant decimation 
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filter, respectively. The up and down arrows represent 
upsampling and downsampling by the associated factor. 
[0055] The input/output relationships of the 

constituent FIR sub-filters ej (n) are given by 

Eq. 10 

N^+Nj-+\-l'\ 
I I 

yM)= Y.^mLAn)^i(n-m) 

where L is the rate change factor, where CmL+i(n) are the 
coefficients obtained from the memory array as described 
above, and where xi(n) and yi(n) are the inputs and outputs 
of the sub-filter e/(n), respectively. In this way, one 
obtains an overall input/output relationship that is 
identical to that of the original cascaded time variant 
and time invariant filters, but the computations occur at 
the lower sampling rate, thus reducing the computational 
cost by a factor of L. Further information relating to 
polyphase decomposition can be found in Vaidyanathan, 
P.P., Multirate Systems and Filter Banks, P T R Prentice- 
Hall Inc., Englewood Cliffs, NJ, 1993, which is 
incorporated herein by reference. 

[0056] As an example, consider a fourth-order combined 
filter with an input/output relationship given by 

Eq. 11 

y(n) = Co (i)x(n) + Cj (i)x(n - 1) + C2 (i)x(n -2) + c^ (i)x(n - 3) + (i)x(n - 4) , 

where i is the index used to address the pre-calculated 
coefficient memory array, where Co(i), Ci(i), Czd), 03(1), 
and C4(i) are the pre-calculated coefficients of the 
combined filter stored in the memory array, and where 
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x(n) and y(n) are the input and output of the combined 
filter, respectively. One possible realization of Eq. 11 
is given in Figure 9. 

[0057] For this example, it is assumed that time 

variant filter precedes the fixed filter. Moreover, it 
assumed that Ty^mATioN = ^ (normally obtained from Eq. 6) . 
Therefore, the combined filter coefficients c^(n)(m e (0,4)) 
are calculated in accordance with Eq. 5 and the 
associated algorithm 1 provided above, and the results 
are stored in a memory array as shown in Figure 10. In 
this way, the corresponding filter outputs can be 
obtained as given in the table of Figure 11. 
[0058] Next, consider the filter in the above example 
when used as an interpolation filter. One possible 
realization of the combined interpolation filter is given 
in Figure 12 . Note that the sample rate of the input 
x(n2) is a factor of three higher than that of the input 
x(ni). As in the example above, it is assumed that the 
time variant filter precedes the fixed filter. Therefore, 
the combined filter coefficients c^(n)(m e (0,4)) are 
calculated in accordance with Eq. 5 and Algorithm 1, and 
the results are stored in a memory array as shown in 
Figure 10. In this way, the corresponding filter outputs 
can be obtained as given in the table in Figure 13 . 
Notably, the computational load is increased 
significantly because each of the additional samples 
added to facilitate the upsampling results in extra 
calculations. However, since the added sample values are 
zero, these calculations are unnecessary. The result is 
a significant, but unnecessary, increase in computational 
load. 

[0059] Finally, consider the use of polyphase 
decomposition of the interpolation filter in the above 
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example as realized in Figure 14. The sample rate of the 
output y(n2) is a factor of three higher than that of the 
input x(ni) . In essence, the number of phases is equal to 
the upsampling factor. The output yfn^J is cycled through 
each phase output for each input x(ni). Accordingly, the 
output y(n2) is effectively upsampled. As in the 
previous examples, it is assumed that the time variant 
filter precedes the fixed filter. Therefore, the combined 
filter coefficients c„ («)(m e (0,4)) are calculated in 
accordance with Eq. 5 and Algorithm 1, and the results 
stored in a memory array as shown in Figure 10. In this 
way, the corresponding filter outputs can be obtained as 
given in the table of Figure 15. 

[0 06 0] Note that the polyphase decomposition of the 
interpolation filter generates outputs that are identical 
to those of the original interpolation filter. Reference 
is made to Figures 13 and 15. However, there has been a 
threefold reduction in computational complexity in the 
polyphase filter, because the adder and multiplier units 
operate at the lower sampling rate and the zero values 
associated with the upsampling are not computed. 
[0061] The following outlines the significant 
reduction in processing provided for the present 
invention as compared with alternative techniques. For 
implementing the regeneration function 74 in a code 
division multiple access (CDMA) time division duplex 
(TDD) system, the following processes are required: 

- spreading the symbol estimates with a real -valued, 
user specific code (length 1-16) ; 

- scrambling the (spread) symbol estimates with a 
complex-valued, cell specific code (length 16) ; and 

- filtering the spread and scrambled symbol 
estimates by a channel estimate (28 - 64 complex taps) . 
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[0062] A graphical representation of a traditional 
regeneration function is provided in Figure 16, wherein 
the symbol estimates are upsam.pled by a spreading factor 
(SF) using upsampling circuitry 76; mixed with the 
appropriate spreading and a scrambling codes using mixers 
78 and 80; and filtered with channel estimates to provide 
a regenerated signal using a complex channel filter 82 . 
Alternatively, the spreading and scrambling process can 
by represented as a multi-rate interpolation filter 84 
with time variant coefficients as illustrated in Figure 
17. Thus, one embodiment of the present invention 
effectively combines the multi-rate interpolation filter 
84 with the channel filter 82 as illustrated in Figure 7. 
[0063] Assume for the purposes of regeneration that 
there are eight users, the spreading code has a spreading 
factor of eight, and channel estimates have a length 
equal to 57. The embodiment of Figure 16 would require 
73 complex words of memory for coefficients and 
spreading/scrambling codes per user and 58 complex 
multiply accumulates (CMACs) per chip-rate period, per 
user, per cancellation iteration. For four cancellation 
iterations, this corresponds to a computational 
requirement of 7.2 billion CMACs per second. 
[0064] In contrast, the present invention, as 
implemented in Figure 7, would create one combined time 
variant polyphase filter for each user, which requires 
1024 complex words of memory for coefficients per user. 
The coefficients are pre-calculated once per slot 
resulting in 0.684 million CMACs per second per user. 
The filter realization requires 8 CMACs per chip-rate 
period, per user, per cancellation iteration for a total 
computational requirement of 0.988 billion CMACs per 
second. The result is functionally equivalent with a 
computational complexity reduction of 86% as compared to 
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the conventional regeneration process. Accordingly, the 
present invention provides functionally equivalent 
results to previous regeneration techniques, while 
significantly reducing the complexity of the chip-rate 
regeneration process. Given that chip-rate regeneration 
accounts for 80-90% of the total computational 
requirement of chip-rate PIC, the invention potentially 
provides more than a threefold reduction in overall 
system requirements for PIC. 

[0065] The invention has many possible applications, 

of which regeneration of signals in TDD MUD is only one 
such application. Those skilled in the art will 
recognize improvements and modifications to the preferred 
embodiments of the present invention. All such 
improvements and modifications are considered within the 
scope of the concepts disclosed herein and the claims 
that follow. 



